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Abstract 

We analyze a second order, linear, elliptic PDE with mixed bound- 
ary conditions. This problem arose as a limiting case of a Markov- 
modulated queueing model for data handling switches in communica- 
tions networks. We use singular perturbation methods to analyze the 
problem. In particular we use the ray method to solve the PDE in the 
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limit where convection dominates diffusion. We show that there are 
both interior and boundary caustics, as well as a cusp point where two 
caustics meet, an internal layer, boundary layers and a corner layer. 
Our analysis leads to approximate formulas for the queue length (or 
buffer content) distribution at the switch. 

Keywords: asymptotics, elliptic PDE, ray method, probability distribu- 
tion. 

MSC-class: 34E20 (Primary) 60J20 (Secondary) 

1 Introduction 

In a model proposed by Anick, Mitra and Sondhi [Tj, a buffer receives mes- 
sages from N statistically independent and identical information sources, 
that asynchronously alternate between exponentially distributed periods in 
the "on" and "off" states. While "on" , a source transmits data at unit rate. 
The buffer depletes through an output channel, with a given maximum rate 
of transmission C. The rate at which a source turns "on" is equal to A and 
the "off" rate is \i. If C < N the buffer may be non-empty, and the condition 



is needed for stability. This simply says that the mean number of "on" 
sources (each transmitting data at unit rate) must be less than the total 
transmission capacity of the channel. This model is analyzed exactly in pQ, 
and the asymptotic limit N — > oo, with 



is studied in [7j . This limit is referred to as "heavy traffic" . 

Analyzing the steady state joint probability distribution of the number of 
active sources and the buffer content, involves solving a system of N linear 
ODEs. In heavy traffic this can be simplified to a backward-forward parabolic 
PDE of the type in (J2J). This model has the disadvantage of treating the buffer 
content as a deterministic fluid. 

A modification of this model, which allows for service variability, is as 
follows. Again there are N independent and identical sources. When a 



A 



N < C 



A + fi 
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source is "on" it generates a Poisson arrival stream to a queue. In the "off" 
state no arrivals are generated. The service time distribution is allowed to 
be general. The model just described may be called a Markov-modulated 
M/G/l queue. 

In jH] it is shown that the joint steady state distribution of the number of 
active sources, the queue length and the elapsed service time of the costumer 
presently being served satisfies a complicated system of integro-differential 
equations. In the heavy traffic limit, where iV — > oo and the average arrival 
rate is close to the mean service rate, this system may be approximated by 
the following BVP: 

Df yy + (c - £)f y + + (£f)s = 0, < ?/ < oo, -oo < £ < oo 

Df y (0,S) + (c- 0/(0, = 0, -oo<£<oo (1) 

oo oo 

J j f(y,Z)dyd£ = l. 

-oo 

Here the variable y is related to the queue length, £ corresponds to a scaled 
measure of the number of "on" sources above their mean value, c > is the 
normalized excess of the service rate over the mean arrival rate, and D > 
measures variability effects in the service time distribution. 

The exact solution to ((TJ) was analyzed in jH] . It is not completely explicit 
and involves finding one eigenvector of an infinite matrix, whose elements 
are complicated expressions involving Laguerre functions. This (infinite!) 
eigenvector must be computed numerically. In the same paper the limit 
D — > oo was considered. Now the matrix becomes diagonally dominant and 
much more explicit results can be obtained. 

The (highly singular) limit D — > was studied in jl], resulting in a 
very complicated asymptotic solution involving contour integrals of parabolic 
cylinder and Airy functions. When D = we see that the problem (JTJ 
degenerates into a parabolic one, that is forward parabolic for £ > c and 
backward parabolic for £ < c : 
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(c - g)Q y + % + (f 9f) { = 0, < y < oo, -oo < £ < oo 

3(0,0 = 0, c<£ (2) 



Now Q is a density in £ and a distribution in y. The problem (J2J) corresponds 
to the heavy traffic limit of the fluid model in pQ. Knessl and Morrison [2j 
derived the exact solution of (J2J). The limit c — ► oo was studied in jH] by 
using the saddle point method and in [9 by using the ray method jS]- 

In this paper we will solve (0) asymptotically in the limit c — > oo by using 
the ray method, the boundary layer method and asymptotic matching j3] . In 
doing so, we shall analyze no less than seven different scales, and one more 
will be briefly discussed in the conclusion section. The asymptotic structure 
of (0) proves much more complicated than that of (0) in the same limit j^]. 

To analyze (0) for large c, it is convenient to introduce the new variables 
77 = £/c, x — y/c, and the small parameter e = c~ 2 . Then JTJ becomes the 
following problem for F(x,rj) = £ _1 /(|/,£) : 

e(DF xx + F m ) + (1 - rj)F x + r]F v + F = 0, x > 0, -00 < rj < 00 

DeF x (0, rj) + (1 - rj)F(0, rj) = 0, -00 < 77 < 00 (3) 



//F(,,,)«fai, = l. 



-00 

The boundary condition together with the normalization condition imply 
that the marginal distribution in 77 is the Gaussian 

00 2 
F{x, V )dx=^ew(-^). (4) 





variable, i.e., 



An important quantity to compute is the marginal distribution in the x 

00 

M (x) =J FMdn . (5) 
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In section 2 we consider the case when x is close to and 77 < 1; this 
will be very useful to match with other asymptotic solutions. Section 3 is 
dedicated to using the ray method to analyze © for e — > with x, r\ fixed. 
This yields asymptotic solutions in two main regions separated by the curve 
x = 77 — hi (77) — 1, 1] > 1. We also derive boundary layer solutions for 

2 

x = 0{£3) and 77 > 1, x — 0(e) and 77 > 1, a corner layer solution in the 
neighborhood of the point (0, 1) and in section 4 a transition layer solution 
along x = 77 — ln(?7) — 1. We show that all the solutions asymptotically match 
to each other in the appropriate limits and also agree with the approximation 
found in section 2. In section 5 we summarize and discuss the main results. 
In section 6 we check the identity (@J) for F(x,r]) and compute the marginal 
distribution in x. 



2 An expansion for small x 

To solve (JHJ) for e small, we will first consider the scaling x = 0(e). Thus we 
introduce the variable v = x/e and convert (J3J) into the problem 

DF VV + (1 - rf)F v + e(r t F ri + F) + e 2 F m = 0, v > 0, -00 < rj < 00 

DF v (0,r]) + (l-r))F(0,ri) = 0, -00 < 7] < 00 (6) 

F(v, 7])dvdr] = -. 

e 

-00 

On this scale (JIJ) transforms to 

J F(v, V )dv=^exp(- 1 £j. (7) 


We consider solutions to © which have the asymptotic form 

F(v }V )~£Le*p(-^\ [F^(v, V ) + ^F^(v, V ) + 0(e)]. (8) 

Substituting (JHJ) into (0) and equating the coefficients of like powers of e we 
get to leading order the equation 

DF^ + (1 - V )F^ = 
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with boundary condition 

DF^(0,T]) + (1-77)^(0,77) = 0, — 00 < 77 < 00. 
Solving for F^°\v,r]) and taking into account (J7|) we conclude that 

1 — 77 



DV2n 



exp 



rf (1 — rj)v 
2~e D 



77 < 1. 



We summarize below the main result of this section. 



Proposition 1 For x = ve 

solution to leading order 



F{v,r}) 



1 — 77 



DV2tt 



exp 



0(e) the equation has the asymptotic 
(1 — r])v 



~2e 



D 



77 < 1. 



(9) 



We see that for 77 = 77^, 77 = 0(1), the solution decouples into a Gaus- 
sian in 77 times an exponential function of v 



F(v,rj) 



'2ir 



exp 



w\ 1 

■ — > x — exp 
2| D 



v 
D 



Such a decoupling was also observed in [Oj, where (jBJ) was analyzed in the 
limit D — > 00 with c fixed. The cases where x is small and 77 > 1 or 77 « 1 
are treated in subsections 3.6 and 3.7. 



3 The ray expansion 

Now we consider solutions of (j2J) which have the asymptotic form 



F(x, 77) ~ e Ul exp 



-\l/(x, 77) 



K(x,rf). 



(10) 



We substitute (fTU|) into Q and equate the coefficients of the lowest power of 
£ to get the eikonal equation for \l/ : 



D {^ x f + [%f + ti{%- * s ) + V x = 0, * t (0, 77) 



77 — 1 
D 
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Equating the coefficients of the next power of e yields the transport equation 
for K : 



DK^ XX + K X + 2DK x m x + + 77^ + 2K V % - rjK x + K = 0, (12) 

3.1 The rays 

We solve (fl"T| by introducing the characteristic curves or rays [x(t), rj(t)], 
written in terms of a parameter t. We first consider rays starting from the r)- 
axis, and impose the initial conditions [x(0), f]{0)} = [0, s]. The characteristic 
ODEs for JTIJ) are: 



dx 

~dt 
drj 

~dl 
dt 

d% 

dt 



-2D^ X +7]-l, x(0) = 

-2^ V -T], T](0)=S 

(13) 



d$> T dx T dri — , \ 2 / T x 2 
dt v dt y 1 y vJ 

From (jUJ) we note that \J/(0,^) = —rf/2, which implies that \1/ [x(0), //((])] = 
^(0,s) = -s 2 /2. 

Setting vp^O, s) = A, ^(0, s) = B and solving (0 yields: 

x = (A - B)e l -(A + B + s)e' t - (2D A + 2A + l)t + 2B + s 

7] = (A - B)e l + (A + B + s)e- t - 2A 
® x = A (14) 
V v = (B - A)e l + A 

m = --{A - Bfe 2t + 2A(A - B)e l - A 2 [D + l)t + AB - ^A 2 + ^B 2 - — . 

The constants A, B can be determined by evaluating the eikonal equation (fTTj) 
at x = (corresponding to t = 0), and also using the boundary condition 
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from (JSJ). This yields 

s — 1 

A = and B = -s or B = 0. 

D 

To decide which value of B is the right one, we take the derivative of \I/ with 
respect to s at t = 

- s = A#(o, s ) = A-f-a?(0, s ) + B^-r](0, s) = B. 
as as as 

Replacing A, B in (JT4*)1 we get: 

x = e t — l-t - - — ^ (s - 1) 

e -t + (D + l)e* - 2. 
r/ = e' + — (s-1) (15) 

1 2t , 2e* - (D + l)e 2t -I. 1 , 

-1 + [4e* - 2(t + 1)](D + 1) - e 2 *(D + l) 2 2 
+ 2D 2 1 ' ' 

For t > and each value of s, the first two equations in ()15|) determine a ray 
in the (x, rj) plane, which starts from (0, s) at t — 0. For s = 1 and s = ^xj, 
we can eliminate t from (|15|) and obtain the explicit expressions 

x = X (rj) = rj — ln(^) — 1, s — 1, 77 > 1 (16) 

1 , . ^ x 1 1 2 
x = — 77 — ln(2 — 77 — -D^), s = — , — < 77 < 



D + l ' v ' D + l' D + l _ ' L> + 1 

For s > -pi-j-, we have both x(t) and r](t) increasing for t > 0. For s = j^f, 
x(t) increases and rj(t) is asymptotic to -prj. 

For s < -pq-j- the rays "turn around" and return to x = for some t* > 0, 
with x(t*) = 0, i](t*) < s. The maximum value in x reached bby the ray 
occurs at t = t rmaY 



t = In 

b x max 111 



-2sL> + D + 2 - 2s + (4s 2 L> - AsD - 8s + 4s 2 + D + 4) 

2(1 - s-L>s) 
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Figure 1: A sketch of the rays in Region I for D = 1. 
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For < s < -^q-j- the ray reaches its maximum in rj at t = t vmax 



D+l 

1 , I - .s 

f = — In 



^7 (^77 max) 2 



1 - s - Ds^ 
1-D 2s + 2D - 2 - s 2 D - sD 2 



sy/(l - s) (1 - s - Ds) 



For s < 0, 7)(t) decreases for < t < t*. 

Solving for s in the ^-equation (fT5|) yields 



e- t + e t -2 + Dr l 
S e-* + (D + l)e* - 2 1 J 

and solving in the x-equation gives 

_ -e f + e~ f + Dt + 2t-Dx 
S ~ -De 1 - e* + e~* + 2Dt + 2i + D ' ^ 

Equating (|T7jl and (fTSJ) we get the implicit equation i? = for the rays, 
where 

R(x, rj, t) = [e~* + (-D + l)e* - 2] a; + (3 - Drj - t - Dt - r))e l + (1 + 1 + rj)e~ 
-A-2t + Dr] + 2tr] + 2Dr]t. (19) 

We sketch several of the rays in Figure H They fill Region I, defined as 

Region I = {x > X = r\ — ln(ry) — 1, r\ > 1} U {x > 0, 77 < 1} . 

3.2 Caustics and cusps 

The Jacobian of the transformation in (JT5J) from Cartesian to ray coordinates 
is 



dx drj dx dr\ 
dt ds ds dt 

= [2{t - 2){s - 1)ZT 2 + {-2t -5s + Us + 2)ZT 1 - s + 2ts + l] e* 
+ [-2{t + 2) (a - 1)ZT 2 + (2t - 2ts + 2 - SsjLT 1 ] e - ' 
+ 8(s-l)Zr 2 + 4(2s-l)ZT 1 
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When J = we can solve for s as a function of t, So = s \j = q 

{-2D - D 2 - 4 + 2Dt + 2t)e 2t + 4(D + 2)e* - 2(2 + D + Dt + t) 

~ (-D 2 - 5D - 4 + 2t + 4Dt + 2tD 2 )e 2t + 8(D + l)e* - 3D - 4 - 2t - 2Dt ' 

(21) 

The equation for the caustic(s), i.e., the points in the (x, 77) plane at which 
the Jacobian is zero, can be given in parametric form. We replace s by Sq in 
the equation of the rays, and let x ca = x(t, S ), r\ ca = rj(t, S ): 

Xca = [—(D + l) 2 e 3t + {2D 2 t 2 - 3tD + D 2 t + 2t 2 - At + D 2 + At 2 D + 6D + 8)e 2t 
-2(3D + 7)e* - e - * + 2(D + l)t 2 + (3D + 4)t + 2(D + 4)] (22) 
/ [{2D 2 t + ADt-A + 2t-D 2 - 5D)e 2t + 8(D + l)e* - (3D + 4) - 2(D + 

^ = [—(D + l) 2 e 3t + 2{2tD + 2t + 2D~ l)e 2t + 2(4 - 2t - 2tD - D)e* + e~ l - 6] 

/ [(2D 2 t + ADt - 4 + 2t - D 2 - 5D)e 2t + 8(D + l)e* - (3D + 4) - 2(D + l)t] 

(23) 

In Figure 121 we sketch the caustic curves for D = 1. There is also a cusp 
where the two caustics meet. Our numerical studies show that the basic 
structure (i.e., the two caustics coming together as a cusp) occurs for all 
D > 0. 

Outside the caustic region, the correspondence between (t, s) and (x, 77) is 
one-to-one. When we are exactly on the caustic curves, the correspondence is 
two-to-one, and inside the region bounded by the two caustics it is three-to- 
one. In Figure El we sketch more densely the rays for D = 1 to indicate this 
correspondence. The evaluation of (fTUj) near caustics and cusps is discussed 
in more detail in section 5. 



3.3 The transport equation 

Now we shall solve the transport equation (|12|) by using (|13|) to write it as 
an ODE along a ray: 

dK 

— = (DV XX + * m + l)K. (24) 



After some algebra, we can show that 

1 1 dJ 

2 ~ 2J ~dt 



D^! xx + * m + 1 = - - — — 



12 



1- 




Figure 3: A sketch of the rays in Region I for D = 1. 
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and hence 



K(x, rf) = k(s) 



62 



To determine k(s) we evaluate the previous result at t — : 

1 



K(0,s) = k(s) 



Using the approximation (JSJ) and the fact that s = r\ at t = 0, we get 



k(s) 



1 — s) 2 , s < 1 and i/i 



3 
2' 



The same result can be obtained by using the BC K x (0,r]) = in (fT^j) 
and fixing the multiplicative constant by normalization. So far we have 
determined \l/ and K only for s < 1. Thus we divide the half-plane x > 
0, — oo < 7] < oo into two parts. The portion filled by the rays for s < 1 
we call Region I and the remainder of the half-plane we call Region II. The 
latter is a shadow of the rays (see also Figure 3.1). 

To summarize, we have established the following. 

Proposition 2 The solution of in Region I is asymptotically given by 



F(x,rj) ~ e 2K(x,i])exp 



-ty(x, rf) 



where 



K(x,r)) 



2,7T 



C2 



V(x,r)) = ~7,e 2t + 



2e l -(D+ l)e 2t - 1 



D 



+ 



-1 + [4e f - 2(t + !)](£> + 1) - e 2t {D + 1) 
2D 2 



(s-lf 



(25) 



(x, rf) is related to (t,s) by / fT3)) and J(t,s) is defined by WDi) . 
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3.4 Region II 

For this region, we consider solutions of (jUJ) which have the asymptotic form 



F(x, rf) ~ e U2 exp 



L(X,T]) 



The term e~^T(x, rj) in the exponent must be included in order for the ex- 
pansion to asymptotically match those valid for small x and r\ > 1, which we 
construct later. 

It follows that $ satisfies (jllj) . L satisfies ()12|) and for T we get the 
following PDE 

(rj - 1 - 2Z>$ !B )r sB - (2$, + V )T V = 0, (26) 

which is equivalent to ^- = 0. Thus we conclude that T is a function of cr 
only and write T(x, rj) = T(a). Here (r, cr) are the new parameters for the 
ray which apply in Region II. Thus a ray starts at r = from rj = o > 1 and 
enters the domain for r > 0. 

The solutions of the characteristic equations are: 

x = (b- a)e T + (a + b- o)e' T + [2a(D + 1) - l]r - 2b + a 
rj = (b — a)e T — (a + b — cr)e~ r + 2a 
= -a (27) 
§ v = (a — b)e T — a 

$ = - a 2 (D + l)r + 2a (a - 6) (e r - 1) - - (a - bf (e 2r - l) + $ (cr). 

Here $o(cr) is the value of $ at r = 0, which corresponds to the ??-axis for 
77 > 1. 

Since from the result for Region I 4| = = for s = 1, we impose 
the condition ^(0, cr) = for all cr > 1. This means that the boundary x = 
will be a caustic curve for rj > 1. Then a has the value 

a(a) = ] —^. (28) 

y ' 2D y ' 

Evaluating (fTTj) at x = we get 

Da 2 + 6 2 + a{b - a) + a. (29) 
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Using (J28)) in (J29)) and solving for b we find that 



o = — ± — — =-, pier = Da + (a — 1) . 



For small r we get from (JTTj) and (}2"%|) 



//9(<r) 



and this implies that the solution b = | — -^yjp must be rejected, in 
that the rays enter the domain x > 0, as r increases. Hence, 



6(a) 



V 7 ^ 



2^ ' 



To find ^ol ") we impose the continuity condition $o(l) = ^(0, 1] 
Since 



d 
da 

we conclude that 



$(0, a) = - a ^-x(0, a) - 6-^-77(0, <r) = -6 



b (u) du 



a- 



4 4 4v^D 



(T 



£ + 1 



+ 



D 



(D + l) 
Dia 



arcsinh 



D 



(D + l)a-l 



D 



arcsinh 



As before, the transport equation (fT2Jl can be solved to obtain 



L(r,a) = L (a) 



(=2 



J 



where 
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J 



dx drj dx drj 
dr da da dr 



-a + l + ~r(a-l) 



D" 



D~2 



+ (-cr -^t + to)d 1 + ^r^/J{a)D-^ + ~ T a\ e T 
+ f-<7 + - raj D" 1 + ^Ty//3(a)D-% - ^raj e~ r 



(34) 



1-1 



+ 2{a - 1)D- Z + 2aD 

In Region II J = only for r = 0. To determine Lo(a) and r(cr) we shall 
analyze the problem for small x, and we will find that not one, but two 
boundary layer expansions are needed to satisfy the boundary conditions (JHJ) 
in this region. 

3.5 Approximation for x = 0(e 5 ) 5 77 > 1 (inner solution) 

We introduce the stretched variable /i = e~ix } and transform (j^J) into 

(1 - 77)^ + + ef ( V F V + F)+ e^F m = 0. (35) 

We represent F in the asymptotic form 

7] — 1 



F ~ e^exp <^ e -1 $ (»7) + e~ 



2D 



-ti + r (77) 



Ro(li,r)) + eaR 1 (ji,n) 



which when inserted into (|33|) give the following PDEs for i?o ; -R 



1 • 



2D 



- ihPt " +[2& ( V )+ V ] [2DT'( V )+^R = 



dfi 2 



= 2D 



2 d 2 R 1 



+ [2V ( V ) + V ] [2DT'(r])+fM} R, 



d/2 2 

2Di[2$M+ V ]^ + [&>( V ) + l]R, 



(36) 



(37) 



(38) 
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Solving (|37j) we get 



R Q = C x {r,) Ai{ 




(39) 



where Ai(-) denotes the Airy function and (3 (rj) is given by (|3*0J). Using (}3T?j) 
into (jHHj) an d solving for Ri we obtain 



Rl = L 2 %D\PW-**CMPW>T? Ai(/Z) + [2L>/3(7 7 )]^ 1 (7 7 )r"(ry)/7Ai(/l)+ 



The function C\{rj) will be determined below. This solution can't satisfy 
the boundary condition (JHJ), and thus we require another boundary layer 
expansion, where x = o(es). 

3.6 Approximation for x = 0(e), 77 > 1 (inner- inner 
solution) 

We introduce the variable v = x/e, and transform (JHJ) to 




with 



71 = 2-*D-*p (77)* [/x + 2^r'(r7)] 



a{n) = (D + 1)77 — 1. 



(41) 



DF W + (1 - 77)^ + efaF, + F) + e 2 F w = 
DF v (0,7i) + (l-ri)F(p,7i) = 0. 



(42) 



We seek solutions of the form 




(43) 



Using (}4*3*j) in (}4*2*j) and taking into account that 
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yields 



dw vv + (i-r))w v + ^(r/) [%{n) + v] w = 

2W t (0, 77) + (1 - ?7)W(0, ry) = 



whose general solution is 



W(v, rj) = w(r]) 



2D 



(77 — + 1 



(44) 



The next step will be finding a corner layer solution valid in a neighbor- 
hood of the point (0, 1), that matches to both the approximation @ and the 
inner-inner solution. This will allow us to determine and w(rj) explicitly 

3.7 Corner layer 

Let us first write F(x, 77) = e V5 exp ( — 3- J G(x, 77), which transforms © into 



(45) 



DsG xx - r]G v + eG m + (1 - rf}G x = 
DeG x (0, V ) + (l-r ] )G(0,r ] ) = 0. 



2 1 

Then we introduce the stretched variables [i = e~^x and 7 = (r/ — 1) 
and (|43j) becomes 



£ 3 G 77 - £ 3 7G 7 + DG W - 7^ - G 7 = 
^ t (0,7)-7G(0, 7 ) = 0. 



(46) 



To leading order G(fx, 7) ~ G(fx, 7) where 



-DG MAt — 7G M — G 7 = 
J DG M (0, 7 )-7G'(0,7) = 0. 



(47) 



The solution to matches to (JHJ) (with /i = 0) if 

^G(0, 7 ) 



1—77 _3 

-e 2 



2ttD 



7 _ 7 

£ 6 , 7 



2ttD 



-00 



(48) 
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so that 1/5 = — ~. In [TO] an explicitly solution to (jjZj) and (f4*Kj) was obtained 
1 , with 



G(/i, 7 ) 



exp 



EL — 7 

2D 12D 



2vr23 J D3 



2vu 



exp{2~3£H 7 A} 



Br 



Ai {\ + 2~hD~l^ 
[Ai(A)] 2 



where Br is a vertical contour in the complex A-plane on which Re(A) > 
and Ai(-) is the Airy function. 

By combining the preceding results we have, on the corner scale, 



F(x, 77) ~ e e exp {^cO, 7)} L c {n, 7) = F(fi, 7) 

T / \ V 2 ^7 7 3 
^ (/i ' 7) = -2l + 2D-12D 



(49) 



1 



exp < 2 3_D 



Br 



Ai (^A + 2-3D-|^ 
[Ai(A)] 2 



dX 



2 1 

where // = e~$x and 7 = e~s (77 — 1) . 

In jlUt Theorem 4] several asymptotic expansions for a function closely 
related to (J4*9*j) were obtained. We use these results in the following sections 
in order to match the different solutions that we have found so far, and 
determine the unknown functions and constants. 



x The function Q(X,T) in [TU] is related to G(ji,j) by 



G(n, 7 ) = ^=(2D)-3Q (2£>)-3 M , (2£>)"5 7 , r = 2* p.. 
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3.8 Matching the solution in Region II and the inner 
solution 

From (|27jl we get the local inversion between (r, a) and (x,r]), for x —* 0, 



r ~ >/2I>^(»7)-^3 + l^p\x + ^(^-Id"? [14/3(77) + HD/Jfa) - 20D] 

o fj{7]) 6b 

(50) 



3 
X - 



a ~ 77 — 

-4 X 2 -f 



1 



/2 

+ ^-/3(r7)~^Z?-i [10/3(7;) + ^(77) - 4£>] 
Using (jSOjl in we obtain 



L ~ L (r})2 4 



^7) 



x 4 + 0(x 4 ), x — > 0. 



(51) 



Expanding for /x — > 00 yields 



i?o ~ C 1 (r ] )/3(r])~24 2-T2D^— yT% 

/ 7T 



x exp 



3 



(52) 



Using (|52|) in (J36|) yields the expansion of F in (}36|) as /i — > 00. By expanding 
77) for small x we see that the exponential parts match automatically 
and the matching of the algebraic factors implies that 



e U2 L ( v )2- 



Hence we have 



D 

W) 



x 4 = e Vz Ci(r])D^ (3(r]) 2 ^2 12 — ^=/z 4 . 

'7T 



^2 = ^3 + 7 
o 

L (r7) = C 1 (r7) J D^/3(r/)^2-i- 



(53) 



7T 
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3.9 Matching the inner and inner-inner solutions 

We take the limit jj, -> in flU and (gUJ) to get 



Ro ~ Ctirj) |Ai 2i J DM(^)^r / (r?) + ^2*D~*P (77) « Ai' 2§D*/3 (77)5 ^(77) /1 

(54) 

and 



D 

6V/%) 



x Ai 

2 

+ 23 

x Ai' 



2 §£> 3/3(77)^^(77) 



+ 



(55) 



1 „ , ,1 



/), ; (//) , ( _ ii)§/3(r7)-*Ci(77) + i^C^a/H 



In order to complete the matching with the inner-inner solution, we must 
have 



e Ui w(v,rj)\ v _ 



From (jUJ), (JHlj) and (JHH|) we conclude that 



^3 + g = ^4 



Ai 



23^)6/3(77)6 r' (77) 



(56) 



C 1 (T7)-2S J D-e/3(7 7 )'Ai' 



23^)6/3(77)6 r' (77) 



10(77) = 23 



D» ) 9(ii)'Ci(r7)--I>*/3(ii)-»C 1 (i 7 ) 



+ - J D6C 1 (77)«/3-6 



Ai' 



2§D5/3(77)5r'(77) 



(58) 
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If we denote by r the smallest (in absolute value) of the roots of Ai, i.e., 
r = max {z : Ai(z) = 0} ~ -2.33810741 
then we have from (|56|) 

T , (Ti) = 2-^D-^P(ri)~K . (59) 
From (JH7J) and ()58j) we obtain an ODE for C\{rj) 

1 a (77) 1 1 a (77) 1 



(7' ( n ) 1 

and a relation between C\{rj) and 10(77) 

1 



C 1 (r/)=0, (60) 



w(77) 

The solution of (jSUjl is 



77—1 



C 1 (77)2fL>i/3(7 7 )sAi , (r ), 77 > 1. 



Ci(7;) = fco(7/-l)^(77H 

with fc a constant to be determined. 



(61) 



(62) 



3.10 Matching the corner and Region I solutions 

From [TIH Theorem 4 (i)] we have the following result valid when /1 and/or 
I7I — > 00 with 7 — y//! — > —00, 



F(fjt,y)~e 6L / (/i,7)exp{^ / (/i,7)} 
*/(^, 7) = -7^75 {t 3 - 18^7 + [7 2 + 6/i] f } - j £ 



27D 

Ll ^ 7) = 3^ if ^ + ^ ~" + 6/i - 2/i 



get 



(63) 
(64) 



We can invert the ray transformation (fTK|) locally when 2 — > 0, 77 — » 1 to 



(65) 



23 



where 



Z = y/Jjl — l) 2 + 6x. 

Using (jnni) in the ray expansion (j2*3j) yields, as (x, rf) — > (0, 1), 



F ~ e 2 K(x, rj) exp 
K 1 1 1 fl 



2-2(77-1)] 



1 1 

2 ' L> 1 27 9 W 



27 



fa " 1) ; 



which agrees with (|U3*|) . 



3.11 Matching the corner and Region II solut' 

From [TIH Theorem 4 (iv) ] we have 

— 7 

7) ~ e~6L n (ii, 7) exp {^//(/i, 7)} 

_29 _1 

2 -2 127/j 4 



-5l 1 



t [Ai'(ro)]' 



7) = — -^? 2 — 7 3 + — M7 — v^/J 

VP ' ' ; 2e ' 12D 2ZT ' 3D ^ 

+ -2*D~3r 7 - 2~s.D _ 5ro v //J 

which is valid when /j and 7 — ► 00, with 7 — ^//I — > 00. 
Combining (fHTjl . (JH3)) and we have 



L ~ /c 2 12 — = 



which agrees with (JSfij) if 



(77 — l)af 



1 3 

r2^2 



[Ai'(ro)r 



.v^dTT 



X 



24 



Since in Region II F(x, rj) ~ e V2 exp e rf) + e 3r(x, 77) L(x,r]),we 

must have 

^2 = (70) 
We use O in (E3), 422) and (PJ) and find that, as (x,rj) -> (0, 1), 



^7) ^ - (^7 - 1) - |(^7 - I) 2 - Y^(^ - X ) 3 + ~ X ) " ^J^ 1 

I» ~ T(l) + -25,0-3^(7? - 1) - 2-sD-lroV^ 
and from (Jl)Hj) we conclude that 



r(i) = o. 



(71) 



We have now determined all the unknown functions from the previous 
sections and these we summarize below 



L(x,r)) = D~i(a-l)-2~i(3{a)"^ 

7T 



q(q) + v^MCDTT) 



r(cr) = 2 - sZr*ro y p rf M 



[Ai'MF 
(72) 

(73) 



i2oOi,77) = (77 -=2-5^(77) 



o(g) + y^yCP + 1) 



2v r D+T 



(74) 



Ai 



x- 



1 _ 5 . 1 



2 aD 6/5(ry)6^ + r 



[Ai'(ro)]' 



'7T 



qfa) + 

d + v/^d+T) 



__VB_ 
2VD+1 



Ai'(r ) 



(75) 
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With (|72|) and (J7HJ) we have completely determined the ray expansion in 
Region II, with ()74|) we have the inner solution (for x = 0(es) and 77 > 1) 
and with (|?5|) we have the inner-inner solution (for x = 0(e) and 77 > 1). 
We have also show that 

4 3 7 

^2 = --, ^3 = -- and v 4 = --. 

4 Transition layer 

Finally we shall find the boundary layer solution near the curve x = X (rj) 
defined by (|TB)l. which separates Regions I and II. We introduce the stretched 
variable u = (x — X )e~3 and (jSJ) becomes 

-27 ? 2 (r ? -l)F w +r ? 2 (^+i 71 )£U/5F^i-[2r ? (r/ - + F u ] e+rfF m e^ = 0. 

(76) 

When s — 1 (a — 1), t — lnfa) (r = lnfa)) and we have 

3 = J [lnfa), 1] = 2 U + In fa) r/+i U - 377 - ~J = 2 J [In fa), 1] = 2j x . 

(77) 



Since 



we should look for solutions of the form 



F ~ ^ exp |-^ 2 - j T( W) 77). (78) 

Using (fTSj) in (|75|l yields for T the equation 

2D 2 j 2 coP(r ] )T w + V 2 D 3 fT v + /3fa) [D 2 f - 2cu 3 ( V - 1)] T = 
whose general solution is 



77) = g ( 7^ J W — exp 



3 



[(277 - 1) (2/V + 2?7 2 - 277 + 1)] j 



Dj J V A? I 27]D 3 j 3 

(79) 

where (7 is a function still unknown. It will be determined in the next section 
by matching with the corner solution. 
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4.1 Matching the corner and transition layer solutions 

Let us first introduce the new variable Q denned by 



(2D)i 



1 

1 



(80) 



From ^UJ Theorem 4 (ii) ] we have the following result, for 7 — > 00, 
Q fixed 



7 7) 



e 6e2eF(/i,7) ~ e 



7 Z6 



where 



(81) 



T d\. 



Ai (^3 A 

The following properties of p(ft) are established in ^U] 
p(0) = 2-i 

p(Q) ~ fitv/7r2~s exp | — ^— j , — > 00 



(82) 



Q2-I 



exp < — 2 srofi?, — > — 00 



[Ai'(r )]' 

In order to match with (|79|) . we first note that 

uj ~ (2D)s(r7 - 77 ^ 1 

thus the right side of (|7S|) behaves as 



(83) 



(2D)"3 



(2D)"3 f) 



J y/2D( V - 1) 

1 



exp 



12D + 1 

8 



ft 3 - ^p?2%D-^{ri-l)e 



: exp 



£6^2^ 

Comparing the above with (|8~Tf we must have 

i/ 6 = -1 



12D + 1 o 1 o 2 „ 1 
8 L> 2 4 r 
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and 



9 



<> :i 1 2D • 1 .. , . . 



6 8 D< 



1 2 

" 

7T 



which implies that 



g{Z) = exp |^ (4D 2 + 6D + 3) j ^2~§p [(2D)s Z 
We conclude by writing the complete transition layer solution in (J78j) 



F~—2-l x ULp 
en V Dj 

u 3 



"2! 


77a; 






exp | 




J 





2£ 2Djei 6 1 ; V^J 



2 ^ 3j 3 [(2^7-1) (2^ + 2^-2,7 + 1)] 

= e^Lxoi^, rj) exp {^x (^, ^7;^)} • 

In the next two subsections we will show that 
solutions in Regions I and II. 



(84) 



matches to both of the 



4.2 Matching the solution in Region I and the transi- 
tion layer solution 

When x is close to X Q , we can invert the equations (fTKJ) to get 



77 

t ~ ln(77), S ~ 1 :(x — X ). 



Using the above in (f23|) we obtain 

V 
2Dj 



1 T 1 

-\[/ rvj — 



2r]D 3 j 



_s _3 77 

£ 2 K ^ S 2 



(z-*o) 3 



DV2irf 

From the definition of Q in ()8U|) we see that 



Q = (2D)-3 



^(c^ + £ 3X ) -^£ 5(77-i; 



(85) 



(86) 
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and thus Q — > ±00 when uo — > ±00. Expanding (J84)) for uj — > 00 and taking 
into account (J82|) we have 

1 M 2 3 - 

This agrees with ([SHjl. since x — X = e^u. 



4.3 Matching the solution in Region II and the tran- 
sition layer solution 

For x — » X we get from (}2Tj) 

T ~ ln(77), cr ~ 1 —(x — X ) 

which when used in fl2*7j) . (J75J) and (J72j) yields 

e 2e' ADeji ! 8r]D 3 jf LV ;/ 

-3 (277 - 1) (2ZV + 2r/ 2 - 277 + 1)] i(x - X ) 3 , 

e~3r~-ie"525£)-sro-?-(z-Xo), (87) 

s-lL ~ -e-^-* i_J_2-¥ f ^ § (x - X ) . 
vr[Ai'(r )] 2 VJi7 

From JHH) and (1821) we find that when — ► — 00 



L x ~ 2"iD~f ( ? j w[Ai'(r )] 2 exp | -2sZr*r ?w 



vj/ l j 



This matches with (IHTjl if we take into account (|77|) . 



5 Summary of results and discussion 

Below we summarize the main results of this section, which consist of the 
asymptotic expansions of F(x, rf) in (jHJ) in the various parts of the (x, rf) 
plane. 
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(A) Region I {x > X = r) - ln(r?) - 1, r) > 1} U {x > 0, 77 < 1} 



F(x,rj) ~ £ 2 i^(x, 77) exp 



-#(a;,77) 



x = e *-l- t -^ + 1 ^-^ +g + e " (a -l), 

n = ^ + e " + ^^- 2 (.-D, 



K(x,r]) 



1 . , 3 ea 
2tt vT 



J = [2(t - 2)(s - 1)£T 2 + (-2* - 5s + 4ts + 2)ZT 1 - s + 2ts + l] e* 
+ [-2(t + 2)(s - 1)£T 2 + (2* - 2ts + 2 - Ssj/T 1 ] e - * 
+ 8(s -1)D- 2 + A(2s -l)D-\ 



T , ^ 1 2t 2e l - (D + l)e 2t - 1 . . 
= --e 2 ' + (s - 1) 

-l + [4e*-2(t+l)](D + l)-e 2 *p + l) 2 
+ 2D 2 1 j • 

2 1 

(B) Corner layer x = fxes , rj — 1 = 753 

7 

F(ar, 77) ~ e~6 exp {* c (a*> 7)} ^c(a*> 7) 



1 1 /• r x ! Ai[A + (2D 2 )"^ 

Lc{»,l) = — r^- / exp (4L>)-s 7 A 1 

V2^(2D 2 )3 2m £ I J [Ai(A)] 2 

C ^' 7J 2e 2D 12D' 
(C) Transition layer x — X = ue^ 



dX, 
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F(x,Tj) ~ e 1 L Xo (uj,r])exp{^x (uj,r];e)} , 



L Xo (x, V ) = VlJ^p 



23 TjUJ 
L>3 J 



j = 2 ( 1 + In (77) 77 + (4 - 3?7 - i 



«"< n > = £i/ 



Br 



Ai 2a A 



,-dA, 



*^" i£) -C-^ 2 + 5( 4B2 + 6B + 3 )(l) : 



2 ^3 j3 [(2r7-l) (2^ + 2^-277 + 1)]. 
Region II {0 < x < X = i] - 111(77) - 1, V > 1} 

4 r 1 1 1 

77) ~ exp 5~ 77) + £~^r(o:, 77) L(x, 77), 



a(a) = — — , b(a) = — + 
V ; 2D ' v ; 2 2^ 

x = (6 - a) e T + (a + 6 - a) e - T + [2a(D + 1) - l]r - 26 + a, 



77 = (6 - a)e T - (a + 6 - (j)e~ T + 2a, 



-a 2 (D + l)r + 2a (a - 6) (e T - 1) - - (a - o) 2 (e 2r - l) + $ 



<7 

$ (a) = ~\- J &(«) du, 
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a(a) = (D + 1) a - 1, /3(a) = Da 2 + (a - l) 2 , 



r = max {2 : Ai(z) = 0} ~ -2.33810741, 



r(tr) = 2 _ 3£rsr y /3(u) * du, 
1 



L(x,?7) = £)-t(<T- l)-2-i/3(a)"^ 

7T 



-D + y/D(D + l) 



2VD+1 



X 



1 ez 



J 



-<7 + l + -r(<7-l) 



D- 2 + 



Iv^)(r-1) 



+ U-T + l)(l-a)D- 2 + 



Iv^)(r + 1) 



+ 2(a-l)D- 2 + 2aD-\ 

2 

(E) Inner layer x = fxes, 77 > 1 



77) ~ £ 2 exp |e 1 



+ e sr (77) ^0(^77), 



V 71 " 



x/TT 
2v r D+T 



Ai 



x 



2 3/} 6/3(77)6^ + r 



[Ai'(ro)] 2 
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(F) Inner-inner layer x = ve, i] > 1 



F(x,r)) 



7 

s exp 



lr) - 1 
2~ D 



-x 



+ e~zT(ri) \ W(v,r)) 



5 1 2 

W(u,7/) = 2~e_£T3 

\/7T 



1 



2VD+1 



Ai'(ro) 



In that part of Region I outside the caustic region (cf. Figure 3.2) the 
mapping between (t, s) and (x, rj) is one-to-one, and K and \I> are unambigu- 
ously determined by the formulas in (A). 

Inside the caustic region the mapping is three-to-one and we should re- 
write (fTUjl as 



_3 
e 2 



A I exp ( -*i j + K 2 exp (-^ 2 j + K 3 exp f-# 3 



where and Kj correspond to the three different values of (t, s) leading to 
the same (x,r]). When t = let us define the starting points on the 77-axis 
of these three rays by the ordering s± < s 2 < S3, where Sj corresponds to tyj 
and Kj. We denote the two caustics by C + and C_ and the cusp where they 

meet clS (^C5 

). Note that the cusp location depends only on D. 

The curve C + has 77 — > —00 as x — > 00, while C_ reaches the ?7-axis 
at some critical point (0, 77*) where again 77* = 77* (D). We have verified 
numerically that along C + we have si = s 2 , \Pi = ^2 an d i^i, -^2 develop 
singularities. However, here ^3 > \&i = ^ 2 and K 3 remains finite. Thus 
we have F ~ e~*K\ exp Q^i) on and near C + . Similarly, along C_ we have 
s 2 = s 3 , \&2 = ^3 and K 2 , K 3 develop singularities. But \&i > \& 2 = ^3 and 

remains finite. Thus the result in (A) remains valid near the caustics, 
except near the cusp point where all three tyj are approximately equal. Here 
the expansion in (A) breaks down. 

Our preliminary results suggest that a new expansion must be constructed 
near the cusp with the scaling 



x - x c = 0(vi), V ~ Vc - A c (x - x c ) = 0(e*). 
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Here A c is the slope at which both C + and C_ hit the cusp. We have thus 
far not been able to complete this analysis. We also note that while the 
expansion near the cusp presents an interesting problem in asymptotics, it is 
not needed for computing the marginal distribution M(x) (0), which is the 
most important quantity from the point of view of applications, and which 
we calculate in the next section. 



6 Marginal distributions 

The last "piece of the puzzle", is to verify that (£Q) is satisfied, and also to 
compute the marginal distribution M(x) in (jSJ). 

We evaluate the integral in (j3J) for e — > 0. For rj < 1, F(x,rj) is concen- 
trated near x — 0, and the result follows from the approximation Q. The 
cases 7] > 1 and rj ~ 1 will be considered below. 

6.1 T] > 1 

In this region F(x,r)) is concentrated near x = X , and using (|51jl and 
we have 



2e 2Dje3 J en V A/ 



2e 2Dje K ' J 2ne V £>J 

and hence, by Laplace's method 



oo 

2 



-oo 

n2 



F(o^)dx~ / ^\-Y £ -2We iX ' Xo) S^\^j dx 



1 / ?y 



exp — — 



This verifies (jlj (at least asymptotically as e — > 0) for rj > 1. 
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6.2 t) « 1 

For ?7 — ^ 1 and x small we use the corner layer expansion, i.e., 
Fix.ri) ~ , \ — rcxp (-— + — - 1 



I /■ r 2 . Ai (a + 2-3,0-3^ 

x / exp ^ 2-3D-37A I ^ '-dX. 



27T« J I J [Ai (A)]' 

Br 

where x = pei and rj — 1 = 755. In the local variable p,® becomes 







so we have to show that 



00 2 
y F(x,77)d/i = £"?-^=exp|-|-| 



A( 7 ) = 2^iexp| I ^| (88) 



where 



A (7) = y ^Lyexp{(A + 2 - to - tAj7 | v 

Br ' 



2 1 



Ai + 



00+i Im(A) 

M h\ I "{^S^ (89) 

Br A 

Taking the derivative of A and using [2] Ai" (p) = p Ai (p) yields 

00+i Im(A) 

(90) 

oo+i Im(A) 



Br A 



1 1 



If f r 21 1 Ai"(p) 



2-3^)3 / / exp i2-3Lr3 7 p} K -^dpdX. 

2m J J P l 7/ 7 [Ai (A)] 

Br A 
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Two integrations by parts give 

oo+i Im(A) 

exp |2-§D _ 57p| Ai" (p) dp 



n 2 d 



i A '< A) l Tx 



exp 1 2 sX) s-yA 

oo+i Im(A) 



Ai(A). 

exp |2~il)~57p| Ai (p) dp 



which when used in (|9U|) leads to the differential equation 

A'(7) = T^7 2 A( 7 ). 



(91) 



Solving (}9T]) yields 



A(7)=Aoexp|^} 



(92) 



where A is a constant. To determine A we let 7 — > —00 in (J89j) . Expanding 
the double integral by a combination of the Laplace and saddle point methods 
leads to 

A( 7 ) ~ D§2^exp . 7^-oo- 

Comparing this to ()92|) we obtain Aq = D%23, which verifies 



6.3 The marginal distribution M{x) 



To evaluate (0) by Laplace's method, we find where \& and $ are maximal 
as functions of r\. We thus examine the equations ^ = and = 0. 

We recall from (|27j) that ^ = (a — b)e T — a. The equation ^ = then 
reads 



VD(<r-l) 



VD(a - 1) + Z}5 ff + D\ I Da 2 + (a - 1) 



< 1, for all D > 0, a > 1 
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We conclude that there is no solution to $^ = for r > 0, and hence < 
in Region II. 

From (JT4"j) ^ v = (A — B)e t — A and consequently 



^ = t = In 



1 - s 



_1 -(£> + !> 



(93) 



which when used in (JTSJ) yields 



^ = <^> ar = Xi(77), < rj < 



1 



£> + l 



Xifa) = -2r]- — (2Dr] - D + 2r] - 2) In 



1 — 77 



1 - (£ + 1)77 



(94) 



The equation x = Xi(rj) defines implicitly r\ as a function of x, 77 = E[x). 
We introduce the function 



^i(x) = ^ [x, E{x)\ 

and from (|15|) we get 

E(x) [1 - £(x)] D + 1 



(95) 



+ 



D ' D 2 

From the defining equation 



[1-E(x)fln 



1- (D + l)E(x) 
1 - E{x) 



-2E(x) + [2(L> + l)E(x) - D - 2] In 

we obtain the asymptotic results 

x 1 x 2 ID - A 



1 - (L>+ l)-B(x) 
1 - E(x) 



x 



E(x) 
E(x) 



x , x — > 



L> 2D 6 D 2 

DTT-(DTi7 exp l- x 



X — >• 00. 



(96) 



(97) 



Use of Laplace' s method to evaluate the integral in (jSJ) as e — >-0 yields 



M(x) ~ e~iK [x, £(x)] V27T 



1 



exp I -\l/ 1 (x) 



vZ-e- 1 ^^,^)] [e 
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and from (|15[) after some algebra we have 



M(x) 



_Al-E(x)f 



exp < -^i(x) 



(98) 



A 



2 [1 - (£> + l)E(x)} [1 - ff(x)] [x + 2ff(x)] (g + 1)D 
2(D + l)E(x) -D-2 
+ D [D + 2£(a;) - 2(D + l)£(a;) 2 ] . 

We can get more explicit results if x is either small or large, using (jHZj). 
We obtain 



M(x] 



D \ D 



M(x) 



x exp 



exp 
D 



1 

e 

+ 



x x 
D + 2LP 



2D + 1 _ r _^ 



x -> 0, 



(99) 



;i + j d) 2 ^(i + d) 2 

x 1 



+ 



1 + .D 



x 



oo. 



The first result in (|9~H|) shows that M(x) is concentrated in the range x = 0(e) 
and the second result is consistent with the spectral solution to (jHJ) obtained 
in 0. 
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